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In magnetic memory and logic devices, a magnet's magnetization is usually nipped with a spin polarized 
current delivering a spin transfer torque (STT). This mode of switching consumes too much energy and 
considerable energy saving can accrue from using a multiferroic nanomagnet switched with a combination of 
STT and mechanical stress generated with a voltage (VGS). The VGS mode consumes less energy than STT, 
but cannot rotate magnetization by more than 90°, so that a combination of the two modes is needed for 
energy-efficient switching. 



In magnetic logic and memory devices, the magnetiza- 
tion of a nanomagnet is usually switched with a spin- 
polarized current delivering a spin transfer torque^. 
This is a well-established technique that has been ex- 
perimentally demonstrated in different systems^ and is 
widely used for spin transfer torque random access mem- 
ory (STTRAM). 

Unfortunately, this method of switching dissipates too 
much energy and a more energy-saving approach is to 
rotate the magnetization of a multiferroic nanomagnet 
- consisting of a piezoelectric layer and a magnetostric- 
tive layer — with an electrostatic potential applied to the 
piezoelectric layer—. The applied voltage generates stress 
in the magnetostrictive layer, which rotates its magneti- 
zation. Such rotations have been experimentally demon- 
strated^. Recently, we showed that this can implement 
Bennett clocking in nanomagnetic logic arrays^ where a 
very small voltage of ~16 mV can rotate the magne- 
tization of a shape-anisotropic multiferroic nanomagnet 
(made of lead zirconium titanate (PZT) and Terfernol- 
D) by ~90° to carry out Bennett clocking. The energy 
dissipated is only ^60 kT at room temperature and the 
rotation takes about 80 ns to complete^. The time can 
be reduced to ~3 ns by increasing the stress to the maxi- 
mum that PZT might be able to generate in the Terfenol- 
D layer, but that also increases the energy dissipation 
to about 48,000 fcT 7 . Still, this is considerably smaller 
than the energy that would have been expended had we 
clocked the same nanomagnet with spin transfer torque 
(STT) in the same 3 ns& 

Unfortunately, the voltage-generated-stress (VGS) 
mode of switching has a major shortcoming. It can ro- 
tate an isolated magnet's magnetization by nearly 90° 
(e.g. from the easy to the hard axis of a shape anisotropic 
magnet) but cannot rotate past 90° to achieve a complete 
flip. Thus, it may be adequate for Bennett clocking which 
requires rotation by ^90°, but inadequate for writing a 
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FIG. 1. An elliptical multiferroic nanomagnet in the y-z 
plane. 



bit in STTRAM, which requires -180° rotation. There- 
fore, we have devised a mixed-mode approach where both 
VGS and STT are employed to rotate a magnet's mag- 
netization by ^180°. We will show that this mixed mode 
switching (MMS) results in considerable energy saving 
at any given switching delay compared to using just the 
STT mode. 

In order to study the MMS paradigm, we have solved 
the appropriate Landau-Lifshitz-Gilbert (LLG) equation 
analytically. Our nanomagnet is an elliptical multifer- 
roic consisting of a 2 nm/3.5 nm thick layer of Terfenol- 
D/nickel (magnetostrictive) and a 40 nm thick layer of 
PZT (piezoelectric) as shown in Fig.[TJ Its major and mi- 
nor axes are 220 nm and 100 nm, respectively, which we 
tacitly assume precludes formation of multiple domains. 
This shape anisotropy causes an energy barrier of 32 kT 
to appear between the two minimum energy states along 
the major axis (easy axis). We assume that the magnet's 
plane is the y-z plane and the easy axis is along the z- 
axis. The angle subtended by the nanomagnet with the 
+z-axis is denoted as 9 (see Fig. 1). 

Solution of the LLG equation yields the time rate of 
change of the angle 9. The derivation can be found in 
the supplementary material, but the final expression is 
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where a is the Gilbert damping constant of the mag- 
netostrictive layer, f2 is the volume of that layer, 
7 = is the gyromagnetic ratio, Bq = 

(/xo/2) M 2 S1 [N d ^yy - N d ^ zz ] (energy barrier due to 
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shape anisotropy), M s is the saturation magnetization of 
the magnetic layer, /i is the permeability of free space, 
Nd-yy and Nd~ Z z are the demagnetization factors in the 
y- and z-directions, respectively, B s 

tress — 

(stress-anisotropy energy) , A s is the magnetostrictive co- 
efficient of the polycrystalline magnetic layer, a is the 
stress, B = Bq + B str ess, s = (h/2e)r)I, I is the in-plane 
spin polarized current inducing STT, n is the spin polar- 
ization of the current, and £ is the angle between the spin 
polarization of the current and the +z-axis. In an actual 
STTRAM configuration, where there will be two magnets 
separated by an insulating layer, the quantity s will be re- 
placed by the quantity s{c(V) — ab(V)}, where b(V) and 
c{V) are dimensionless voltage-dependent quantities 8 , 
but for an isolated magnet, they do not arise. 

In order to flip the magnetization from -z to the +z 
direction, the spins in the spin polarized current must be 
aligned along the +z-direction so that £ = 0. 

From Equation ([1} , we can obtain the time r required 
to rotate the magnetization from an initial orientation 9\ 
to a final orientation 62 as 



t= / dt = 
Jo 



1 



8l 7 [s sin{9) + aBsin{29)] 
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We will assume that the initial orientation is aligned close 
to the -z-axis so that dinitial = 180° — e. If e = and 
the magnetization is exactly along the easy axis, then no 
amount of stress or spin polarized current can budge it 
since the effective torque exerted on the magnetization by 
either stress or spin polarized current will be exactly zero. 
Such locations are called "stagnation points" . Therefore, 
we will assume that e = 1°. This is not an unreasonable 
assumption since thermal fluctuations can dislodge the 
magnetization from the easy axis and make e — > 1°. 

We now consider three different scenarios: pure stress 
mediated switching (VGS), pure spin transfer torque 
mediated switching (STT) and mixed mode switching 
(MMS). 

In the pure VGS switching mode, the time taken to 
rotate the magnetization from the initial orientation 9 = 
180° - e to 9 = 180° - (</> < 90°) is obtained from 
Equation ((2|) as 
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The above equation shows that stress will take infi- 
nite time to complete the rotation if either e = 0, or 
<f> = 90°. Thus, initial alignment along the easy axis or 
final alignment along the hard axis is forbidden. The 
physics behind the stagnation points can be understood 
by looking at the energy profile of the nanomagnet E{6) 
as a function of 9 (see the supplementary material for an 
expression that relates E{ff) to 9). Without spin transfer 
torque, the energy consists of only the shape anisotropy 



energy and the stress anisotropy energy. The slope of 
this energy profile E{9) at 9 = 0°, 90° and 180° is exactly 
zero, which means that there is no torque acting on the 
nanomagnet at these orientations to cause any rotation. 
This ensures that the pure VGS mode can never rotate 
the magnetization past 90°. It can rotate the magnetiza- 
tion from close to the easy axis to close to the hard axis, 
which is adequate for Bennett clocking of nanomagnetic 
logic£, but not adequate for writing bits in MRAM. 

Next, let us consider the pure STT switching mode 
with no stress present. The time required to rotate the 
magnetization from an initial orientation 9 — 180° — 1? to 
a final orientation 9 = d is obtained from Equation (|2|) 
as (see supplementary material) 
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where m = 2aB^/s and "d is an arbitrary angle. The 
above expression is invalid for rn > 1 since then the en- 
ergy imparted by the spin transfer torque will not be 
sufficient to overcome the energy barrier due to shape 
anisotropy and cause rotation of the magnetization. 

The last expression shows that the time to complete 
the rotation approaches infinity if $ — > 0, indicating that 
spin transfer torque has a stagnation point at 9 = 0° 
and 180°, i.e. when the final or the initial state of the 
nanomagnet is along the easy axis. However, there is no 
stagnation point at the hard axis (9 — 90°), unlike in the 
case of VGS. Hence, thermal fluctuations are needed to 
nudge the magnetization slightly away from the easy axis 
so that spin transfer torque can nearly complete the flip. 
Once again, the origin of the stagnation points can be 
understood by noting that no spin transfer torque acts 
on the nanomagnet at these locations. 

If thermal fluctuations deflect the nanomagnet from 
the easy axis by an angle e, then the time At needed 
to complete a near-flip is obtained from Equation §Q 
by setting d — e. Note that At becomes exponentially 
longer as e approaches zero, so that one can nearly, but 
not completely, flip magnetization with STT. 

For the MMS mode, we consider the following switch- 
ing sequence which results in the most energy saving. 
Stress is applied abruptly to initiate the switching and 
rotate the magnetization from its initial orientation along 
180° - e to 9 = 180° - $ (■& < 90°), then a constant spin 
polarized current is turned on instantaneously. There- 
after, the sign of the stress is reversed abruptly as soon 
as the hard axis is crossed, and then the spin polarized 
current is turned off instantaneously when 9 — •&. Fi- 
nally stress is removed when 9 = e. This completes the 
switching by rotating the nanomagnet from 9 = 180° — e 
to 9 = e, where e = 1° in our simulations. From Equa- 
tion @, the time required to rotate from 180° — z? to •& is 
tMMS = tmms+ + I'm ms-, where tmms+ and tmms- 
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are given by 



TMMS± = ± 
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-2m±ln\l ± m±cos'd\ 



l{m 2 ± - 1) 

+ (m± - l)^n|l=F cos 0\ + (m± + 1) ln\l ± cos (5) 

where 77i± = 2a (£? ± B stress ) /s. The upper sign de- 
notes the switching delay for the rotation from 180° — fl 
to 90° and the lower sign denotes the switching delay for 
the rotation from 90° to ■&. We add to this the time to ro- 
tate with stress alone between 9 = 180°-e to = 180°-tf 
and from 9 = d to 9 = e which are obtained from Equa- 
tion flS). 

Our next task is to compute the energy dissipations in 
the different switching modes. Since the pure VGS mode 
cannot rotate past the hard axis and flip the magneti- 
zation, we shall not consider it further and concentrate 
instead on the pure STT mode and the MMS mode. 

In the pure STT mode, we will first pick a switching 
delay At and then back out the current I (At) needed to 
switch in that time from Equation ([?]) using appropriate 
material values for nickel and Terfenol-D given in the 
supplementary material. The energy dissipated is then 
computed as Edi SS (At) — I 2 (At)RAt, where R is the 
resistance of the magnetostrictive layer. This allows us 
to plot Ediss versus At for the pure STT mode. 

In the MMS mode, the energy dissipated is 
E Aiss (a,At,ti) = (5/2)CV 2 (a) + I 2 (At,a,tf)RAt, where 
C is the capacitance of the piezoelectric layer (estimated 
as 3.8 fF for the dimensions of the multiferroic nanomag- 
net and assuming a relative dielectric constant of 1000 
for the PZT), V(a) is the voltage applied on the multi- 
ferroic to generate a given stress a, and I(At, a, $) is the 
spin polarized current needed to rotate the magnetization 
from 180° — $ to $ and complete the switching in time 
At in the company of that given stress. We pick a stress 
a and a At, and then back out the current I (At, cr, 
from Equation ([5]) that will be needed to complete the 
switching in time At while applying STT in the angle 
interval [180° -&,&]. Clearly, E diss (a, At, d) depends 
on d. Therefore, we have to first find E,a ss (a, At, •&) for 
given At and cr as a function of •& (see supplementary 
material for a plot of Edi ss (c> At, $) versus •&). We then 
choose that value of $ for which the E,a ss (a, At, $) is 
minimum since that will yield the most energy efficient 
switching scheme. For this optimum value of we find 
Edi as (o~, At) and repeat this procedure for every a, At. 

In Fig. [3J we plot the energy dissipations for the pure 
STT mode and the MMS mode as a function of delay in 
both Terfenol-D/PZT and nickel/PZT multiferroics with 
stress as a parameter. Typical switching delays of mag- 
nets contemplated for STTRAM is 2-20 ns- , which is why 
we plot within this range of delay. In this range, MMS is 
always superior to STT in the same material, and the en- 
ergy saving increases with increasing delay since a larger 
portion of the rotation role can be handed over to stress 
if we allow for a larger delay. The saving is very pro- 
nounced in Terfenol-D/PZT (3 orders of magnitude at 



long delays) but much less pronounced for nickel (36% 
at long delays of 20 ns). However, for delays exceeding 
6 ns, the MMS mode in Terfenol-D/PZT is more energy 
efficient than both the pure STT in nickel and the MMS 
mode in nickel/PZT. 
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FIG. 2. Energy dissipation as a function of switching delay At 
for various applied stress a in the MMS mode for (a) Terfenol- 
D/PZT and (b) nickel/PZT. The plot for the pure STT mode 
is also shown. 

In conclusion, we have shown MMS mode switch- 
ing can be considerably more energy-efficient than STT 
mode. This could potentially lead to a new family of 
magnetic logic and memory devices based on multifer- 
roics that are switched with a combination of voltage- 
generated stress and spin polarized current. 
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SI Magnetization dynamics of a free two-dimensional nanomag- 
net: Solution of the Landau-Lifshitz- Gilbert equation 

Consider a two-dimensional isolated nanomagnet of ellipsoidal shape lying in the y-z plane with 
its major axis aligned along the z-direction and minor axis along the y-direction (Fig. 1 of the main 
paper). The dimension of the major axis is a and that of the minor axis is b, while the thickness 
is I. The volume of the nanomagnet is 0, = (n/tyabl. The magnetization of the magnet always 
lies in the y-z plane since the magnet is two-dimensional. Let 9(t) be the angle subtended by the 
magnetization with the +z-axis at any instant of time t. 

The total energy of the single-domain nanomagnet is the sum of the uniaxial shape anisotropy 
energy and the stress anisotropy energy: 

E = Esha + Est a, (SI) 



where Esha is the uniaxial shape anisotropy energy and Esta is the stress anisotropy energy. The 
former is given by 

Esha = (^ /2)M^N d , (S2) 
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where M s is the saturation magnetization and N d is the demagnetization factor expressed as 

N d = N d _ zz cos 2 0(t) + N d _ yy sin 2 9(t) (S3) 



with N d -zz and N d _ yy being the components of N d along the z-axis and y-axis, respectively. When 
I < a, b, N d _ zz and N d _ yy are given by [S3] 

^ = I (a) 

Note that uniaxial shape anisotropy will favor lining up the magnetization along the major axis 
(z-axis) by minimizing Esha, which is why we will call the major axis the "easy axis" and the minor 
axis (y-axis) the "hard axis". We will assume that a force along the z-axis (easy axis) generates 
stress in the magnet. In that case, the stress anisotropy energy is given by 

Est A = -{3/2)\ s ancos 2 6{t), (S5) 

where (3/2) X s is the magnetostriction coefficient of the magnet and a is the stress generated in it 
by an external agent. Note that a positive \ s o product will favor alignment of the magnetization 
along the major axis (z-axis), while a negative \ s a product will favor alignment along the minor axis 
(y-axis), because that will minimize Est A- In our convention, a compressive stress is negative and 
tensile stress is positive. Therefore, in a material like Terfenol-D that has positive A s , a compressive 
stress will favor alignment along the minor axis, and tensile along the major axis. The situation 
will be opposite with nickel that has negative A s . 

At any instant of time, the total energy of the nanomagnet can be expressed as 

E{t) = Bsin 2 9(t) + C (S6) 
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f^M 2 i 
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Bstress = (3/2)\ s an (S7b) 



B = y Mfr [N d _ yy - N d _ zz ] (S7a) 
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B = -Bo + B stress (S7c) 

C = ^-M^flN d . zz - (3/2)X s an. (S7d) 

Note that Bq is always positive, but B stre ss can be negative or positive according to the sign of the 
\ s a product. 

The magnetization M(t) of the magnet has a constant magnitude at any given temperature but 
a variable direction, so that we can represent it by the vector of unit norm n m (i) = M(i)/|M| = e r 
where e r is the unit vector in the radial direction in spherical coordinate system represented by 
(r,6,4>). The other two unit vectors in the spherical coordinate system are denoted by e# and for 
9 and 4> rotations, respectively. 

The torque acting on the magnetization within unit volume due to shape and stress anisotropy 

is 

T E (t) = -n m (t) x VE[6(t)] = -2Bsin6{t)cos6(t) e (S8) 

since VE[0(t)] = (dE(t)/dO(t))e e and dE{t)/d6{t) = 2BsinQ(t)cos6{t). This immediately shows 
that the torque has out-of-plane component. 

Passage of a constant spin-polarized current / though the magnet generates an additional spin- 
transfer-torque given by [921 92] 

T STT (t) = -s n m (t) x [an m (t) + bn s (t) + e(n m (i) x n s (t))} (S9) 

where s = (h/2e)rjl is the spin angular momentum deposition per unit time and r] is the degree of 
spin-polarization in the current /. In order to minimize the resistance in the path of the current /, it 
will be applied in-plane [91] rather than perpendicular-to-plane, which has been the tradition [95] • 
The coeffcients b and c are voltage-dependent dimensionless terms while a is somewhat irrelevant in 
this context since the term involving a vanishes. The unit vector n s is in the direction of the initial 
spin polarization of the incident current and lies in the y-z plane. For the sake of simplicity, we will 
assume that it is time-invariant. If the current / flows in the y-z plane and the spin polarization 
subtends an angle £ with the positive z-axis, then the spin-transfer torque is given by 

T S TT(t) = s l-b(V) sin(i - 6{t)) e> + c(V) sin(£ - 0(t)) eg] (S10) 

where b(V) and c{V) are the voltage-dependent coefficients [939Z]- 

The magnetization dynamics of the single-domain magnet under the action of various torques 
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is described by the Landau-Lifshitz-Gilbert (LLG) equation as follows. 



dn m (t) / dn m (t) \ 7 

_ + a n m(( )x_ =_ (TE(t) + TsTT(t)) (S11 ) 



where a is the dimensionless phenomenological Gilbert damping constant, 7 = 2/x#//i is the gyro- 
magnetic ratio for electrons and is given by 2.21 x 10 5 (rad.m).(A.s) -1 , and My = [i§M s Vl. In the 
spherical coordinate system, 

dn m (t) 

= e'e e + sin9(t)(j)'e^. (S12) 



dt 

where ()' denotes d()/dt. Accordingly, 

^ in 

n^(t) x — 

dt 



a 



dn m (t) 

n m (t) x — = -asin6(t) 4>'{t) e e + aO'{t) e (S13) 
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and 



dn m (t) 



dn m (t) 

n m (t) x — I = (0'(t) - asin9{t) </>'(£)) e e + (sin9(t) 4>'{t) + a9'(t)) e^. (S14) 
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Equating the eg and components in both sides of Equation (jSlip . we get 

0'(t)-asinO(t)<t>'(t) = -^-s c(V) sin(£ - 0(t)) 

My 

sinO(t)<l>'(t) + a6'(t) = (-2Bsin9(t)cos9(t) - sb(V) sin{£ - 6{t))) . 

My 

(S15) 

Eliminating (j>'(t) from the last two equations, we get the following equation involving only 6(t): 



(1 + a 2 ) 9\t) = -jL. [s c(V) sin(£ - 6(t)) + a {-2Bsin6{t)cos9(t) - s b(V) sin(£ - 6{t)))} 

<s {c(V) - a b(V)} sm(£ - 9{t)) - 2aBsin9{t)cos9{t)] . (S16) 



My 



In the rest of the derivation, we will assume that the initial magnetization of the nanomagnet 
is close to the negative z-axis, so that the initial value of is close 180°. Our intent is to switch 
it, i.e. flip it up, so that the final value of 9 will be close to 0°. Accordingly, the spin polarization 
of the current needs to be along the positive- z direction, which means that £ = 0°. Note that we 
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said that the initial and the final orientations of the magnetization are "close to" but not "exactly 
along" the easy axis, since the torques Te and T$tt both tend to zero when 9{t) tends to either 
or 7r as can be easily verified from Equations ([S8]) and ()S10p . Since the torques become vanishingly 
small near the easy axis, neither the initial nor the final orientations can be exactly along the easy 
axis. In fact, if the initial orientation is exactly along the easy axis, then no torque acts on the 
magnetization vector and nothing can budge it from this location. The locations where the torque 
vanishes are called "stagnation points" . Clearly, the easy axis is a stagnation point for both VGS 
and STT mode switching. Note that the torque Te is zero when 6(t) = it/ 2 as well, so that the 
hard axis is also a stagnation point for the VGS mode, but not the STT mode. We will revisit all 
this later. 

a Switching delay for stress-mediated rotation or for the pure VGS mode 



In the absence of spin-transfer torque, Equation (|S16p becomes 



(1 + a 2 ) e'(t) = J— 2aBsin6(t)cos9(t). (S17) 

which yields 

J dt ~~ J 2 ia Bsine{t)cose{t) (S18) 

Integrating, we get that the time t taken for the magnetization angle 8(t) to change from Q\ to 
02 under the action of stress is 



(1 + a 2 ) n M s Q 

Wgs = 7, o ln 

27 an 



tanOo 



tanO\ 



(S19) 



The above expression shows immediately that we can approach, but never actually reach, the 
hard axis with just stress acting on the magnet (i.e. in the pure VGS mode) since tycs ~^ 00 
as 02 — > 7r/2. In other words, it will take infinite time to reach the hard axis starting from any 
orientation. It should be also apparent that it will take infinite time to get out of alignment with 
the easy axis if the initial orientation was along the easy axis since tycs — >• 00 as — > 7r. Therefore, 
clearly, both easy and hard axes are stagnation points. The switching time tends to infinity in these 
cases since the torque vanishes at the stagnation points, as we have seen already. As a result, the 
VGS mode may be good for Bennett clocking [S8J which merely requires rotating magnetization 
from near the easy axis to near the hard axis, but it is useless for flipping magnetization needed in 
spin valves and magnetic tunnel junctions, since it is impossible to rotate past the hard axis and 
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complete a near-180° rotatio 



b Switching delay for spin-transfer torque (STT) 

In absence of stress, Equation (|S16p becomes (assuming £ = 0) 



'l + a 2 )6'(t) 
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[s' sin9(t) + 2aB sin9(t)cos6(t)] 



(S20) 



where s' = s{c(V) — ab(V)} and B = Bq since B streS s = 0. The last equation can be integrated 
using the identity 



d0 



1 



sin9(t) + msin9{t)cos6 2(m 2 — 1) 

x [(m - l)Zn(l - cos9(t)) + {m + l)/n(l + cos9(t)) - 2mln\l + mcos9(t)\] . (S21) 



which yields the switching delay associated with rotating the magnetization from 9 = TT — -dto9 = r d 
as 

(1 + a 2 ) fi M s Q m 
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(S22) 



where m = 2cxBq/s' . The above equation is valid only for m < 1. If m > 1, then the spin transfer 
torque energy is insufficient to overcome the energy barrier due to shape anisotropy and cause 
rotation past the hard axis. 

Figure [Si] plots the time it takes to rotate the magnetization from 9 = n—fl to 9 = •& as a function 
of the angle i9. It is obvious from Equation ()S22|) that the delay should increase rapidly as the initial 
oientation of the magnetization approaches the easy axis since t$TT diverges when •& approaches 
or 7r (because of the last logarithmic term in Equation (|S22p . This feature is apparent in the plot 
as well. Once again this happens because the spin transfer torque vanishes when the magnetization 
is exactly along the easy axis. Therefore, the easy axis presents a stagnation point and thermal 
fluctuations will be needed to get out of it. However, there is no stagnation point around the hard 
axis (unlike in the case of VGS mode switching) which can be verified from Equation (|S20j) as well. 
This is what allows STT mode switching to rotate the magntization past the hard axis. 

The reader can easily ascertain that the time for flipping from negative to positive z-axis is the 
same as that for flipping from the positive to the negative z-axis. In the latter case, we merely have 
to make £ = it instead of 0, and interchange it — 1? and $. 



1 Bennett clocking in nanomagnetic logic arrays is also greatly facilitated by the dipole interaction of near neighbors 
which comes to the aid of stress. Discussion of this issue is however outside the scope of this paper. 
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c Switching delay for the mixed mode 

The following considerations come into play in choosing the mixed mode of switching. The 
VGS mode cannot rotate past the hard axis and therefore will need assistance from the STT 
mode to accomplish this feat. Therefore, the spin polarized current should be turned on when the 
magnetization approaches the hard axis and turned off after crossing the hard axis, i.e. the current 
/ should be on between 9 = tt — •& and 9 = But how should one determine the optimum #? 
Clearly there is an optimum value of i? for a given switching delay since if $ is too large and close to 
90° , then we would have burnt up a lot of the time with stress which is slow and now we will need 
to pass a very high spin polarized current / to complete the switching in the little time remaining. 
This will be counter-productive since the large current will dissipate too much energy. On the other 
hand, if i? is too small, then we are not being sufficiently frugal in the use of energy-inefficient spin 
transfer torque and once again expending too much energy because we are running the current / 
for too long. This too is counter productive. Hence, there is an optimum value of $ for a given 
delay that must be found from exhaustive search. If the delay is too short, then $ may actually 
turn out of be zero, or rather e, i.e. the spin polarized current must be on for the entire duration to 
achieve the switching within the stipulated delay. A plot of i? versus the switching delay At is given 
in Fig. IS2I Note that for short delays, # is indeed e ( = 1°), which means that the spin polarized 
current is on for the entire switching duration, and then •& increases sharply with increasing delay. 
Clearly, the total energy dissipated in the switching process will depend on •& since that determines 
the duration for passing a spin polarized current. For a given total delay At and a given stress a, 
one must compute the total energy dissipated in both the spin polarized current and the stress as 
a function of i? and choose that value of i? for which the total energy dissipated is minimum. 

The stress too must be handled judiciously in the mixed mode. One should remember that a 
positive X s a product favors aligning the magnetization along the easy axis while a negative product 
favors alignment along the hard axis in the configuration we choose. Since we are starting off from 
a location close to the easy axis towards the hard axis, we should start off with a negative product 
(compressive stress for a material with positive magnetostriction and tensile stress for a material 
with negative magnetostriction) and then reverse the stress to make the product positive once the 
rotation goes past the hard axis, because then we will be aiming for the easy axis. This will be the 
optimum switching strategy. 

Switching via the MMS mode is therefore accomplished in four distinct phases: in phase 1, stress 
of appropriate sign (compressive or tensile) is applied to rotate from 9 = tt — e to 9 = tt — then 
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in phase 2, STT is added to stress to rotate from 8 = tt — $ to 9 = n/2, then in phase 3, the sign 
of the stress is reversed while still maintaining STT to rotate from 8 = tt /2 to 8 = i9, and finally in 
phase 4, STT is turned off and stress alone rotates the magnetization from 9 = •& to 9 = e. This 
nearly completes a flip. 

The time taken to complete the first phase of the rotation from 9 = tt — e to 8 = tt — # is found 
from Equation fjS 19|) 

h 4- cv 2 ) HnM.O i.mn fi 

(S23) 



(l + a 2 )fi M s fl r 
27 a [B + B stress \ 



tane 



Note that in the above equation, # > e. Therefore, to ensure that t v hase\ is positive and hence 
meaningful, we have to ensure that B s t ress is negative (i.e. the A s cr product is negative) and that 
I B s t ress I > -Bo- The former condition determines the sign of stress (compressive or tensile) and the 
latter condition ensures that the stress anisotropy energy can overcome the shape anistropy energy 
and initiate rotation. 

The time taken to complete the final phase of the rotation from 9 = 1? to 9 = e is also found 
from Equation (15191) : 

I i 4- n 2 ) n,,M O tane 



t (1 + a 2 ) noM s n 

1 27 Q [B Q + Bstress] 



tand 

Clearly t p h ase 4 will become negative (unphysical) if B s tress is negative and exceeds Bq in magni- 
tude. To prevent this, the sign of the stress must be opposite to that during phase 1 so that B s tress 
positive. That will keep t p h ase 4 positive. Interestingly, a negative B s tress is still physical (and keeps 
t p hase4 positive) as long as its magnitude is less than that of Bq. In that case, the shape anisotropy 
will ensure that the magnet comes to rest along the easy axis in spite of opposition from the stress 
anisotropy energy because shape anisotropy energy is stronger. Nonetheless, the opposition from 
the stress anisotropy energy will slow the rotation. In order to speed up the rotation, it is therefore 
necessary to reverse the stress during phase 4 so that shape anisotropy and stress anisotropy energy 
work in concert instead of opposing each other. 

During the second and third phases of the rotation, both VGS and STT modes coexist. The 
only difference between these two phases is that during the second phase, B stress = — \B s tress\, while 
during the third phase, B stress = \B s tress\- 

In the presence of both stress and spin-transfer torque, Equation ()S16[) becomes (assuming £ = 0) 

(1 + a 2 ) 9' = J— \s' sin9 + 2aBsin9cos9] (S25) 

HoM s il 
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or 

dt = - [ (1 + M ° Ms ° I ( 1 + q2 )^ M ^ de (S26) 

7 7 [s' sinO + 2aBsin6cos9] J 7 [s' sinQ + <xBsm(20)] 

where s' = s {c(V) — a b(V)}. This is Equation (2) in the main paper. 

The above equation immediately shows that stress (contained in the second term in the denom- 
inator) is most effective when 6 = 45°, while the spin transfer torque (contained in the first term in 
the denominator) is most effective when 6 = 90°. The stress term actually vanishes when 9 = 90° 
consistent with the fact that stress is ineffective around the hard axis. 

We will carry out the integration using the identity in Equation ()S21j) and arrive at the following 
expression for the switching delay: 

tphase2 = T MMS+ 

tphase3 = TMMS— (S27) 



with 

(l + a 2 )n M s n m± 



T MMS± 



2-ya(B =F \B stress \) m\ - 1 
[(m± — 1) ln\l =F cos i?| + (m± + 1) ln\l ± cos $\ — 2m±ln\l ± m±cos (S28) 



and m± = 2a(B =f \B stress \) / s' . 

S2 Material parameters 

The material parameters that are used in the simulation are given in the Tables [T1 and 121 [S 91-9 1 3 j . 



Table 1: Material parameters for nickel 



Young's modulus (Y) 


2xlO n Pa 


Magnetostrictive coefficient ((3/2)A s ) 


-3xl0" 5 


Saturation magnetization (M s ) 


4.84xl0 5 A/m 


Gilbert damping (a) 


0.045 


Resistivity (p) 


7.8xl0~ 8 ohm-m 



Table 2: Material parameters for Terfenol-D 



Young's modulus (Y) 


8xl0 lu Pa 


Magnetostrictive coefficient ((3/2)A s ) 


+90xl0~ 5 


Saturation magnetization (M s ) 


8xl0 5 A/m 


Gilbert damping (a) 


0.1 


Resistivity (p) 


63.1xl0~ 8 ohm-m 
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53 Procedure for determining the voltage required to generate a 
given stress in a magnet ostrictive material 

In order to generate a stress a in a magnetostrictive layer, the strain in that material must be 
e = cr/Y, where Y is the Young's modulus of the material. We will assume that a voltage applied 
to the PZT layer strains it and since the PZT layer is much thicker than the magnetostrictive layer, 
all the strain generated in the PZT layer is transferred completely to the magnetostrictive layer. 
Therefore, the strain in the PZT layer must also be e. The voltage needed to generate this strain is 
calculated from the piezoelectric coefficient of PZT. 

We assume that the maximum strain that can be produced in the PZT layer and transferred 
to the magnetostrictive layer is 500 ppm. This gives the maximum stress in Terfenol-D as 40 MPa 
and in nickel as 100 MPa. 

54 Simulation results 

Some additional simulation results and corresponding discussions are given in the Figures IS3I 

-ES 
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Switching delay for spin-transfer torque 




a (in degree) 

Figure SI: Plot of the time taken by STT to rotate magnetization through an angle $ versus '&. 
The material is Terfenol-D. The initial orientation is assumed to be along the easy axis. Note that 
slope is high near i? = 0, indicating that it takes a very long time to rotate even by a small angle 
away from the easy axis since that location is a stagnation point. 
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STT active during (71-a) to a for stress 40 MPa 




Delay (ns) 

Figure S2: Plot of the angle i? as a function of the total switching delay At in the MMS mode 
for Terfenol-D for a stress of 40 MPa. The angle $ asymptotically approaches 90° for very long 
delays, but can never actually reach 90°, because then no STT mode would have been required to 
accomplish the switching. That is impossible since the VGS mode by itself cannot rotate past the 
hard axis and complete switching. 
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Energy Dissipation with STT+Stress 




Delay (ns) stress ( MPa ) 

Figure S3: Energy dissipation at different delays and at different values of applied stress for the MMS 
mode, assuming that the material is Terfenol-D. These quantities are calculated for the optimum 
angle $ to get the minimum energy dissipation. 
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Energy dissipation with 40 MPa stress and delay S ns 

.65 i 1 1 1 1 1 1 1 1 



0.6 - 



h- 0.55 - 
2 0.5- 




10 20 30 40 50 60 70 80 90 

& (degree) 

Figure S4: Energy dissipation for a fixed delay of 8 ns and for a fixed stress of 40 MPa in Terfenol-D 
as a function of the angle # in the MMS switching mode. 
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Figure S5: STT current required to switch using the MMS mode as a function of switching delay 
for different applied stress in the case of Terfenol-D. The current required to switch using the pure 
STT mode is also shown for both Terfenol-D and nickel as a function of delay. Note that the latter 
current is almost the same for both Terfenol-D and nickel at any delay within the range shown. At 
any given delay, the current needed to switch using the MMS mode decreases with increasing stress 
because stress carries out part of rotation and relegates less and less of the rotation task to STT. 
At high enough stresses there are sudden jumps in the current when the delay exceeds a certain 
value. At these delay values, the angle $ is no longer stuck at e and jumps to a large value (see 
Fig. IS2p so that STT need not be kept on during the entire duration of rotation. It becomes more 
economical to turn on the STT current only for part of the rotation (from 9 = 180° — $ to 9 = •&), 
but that also requires a larger STT current, which is why we observe the jump. 
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Energy dissipation for delay S ns 




25 30 35 40 

Stress (MPa) 



Figure S6: Total energy dissipated (in stress and spin polarized current) to switch the magnetization 
of a Terfenol-D nanomagnet using the MMS mode in a fixed delay of 8 ns as a function of applied 
stress. As we increase stress, we require progressively less spin transfer current to complete the 
switching in the stipulated 8 ns. This results in overall energy saving since spin polarized current 
is far more dissipative than stress. Energy is saved by maximizing the proportion of stress and 
minimizing the proportion of spin transfer torque current when the switching delay is fixed at 8 ns. 



17 



Roy, K. et. al, APL (2010) 



Energy Dissipation with STT+40MPa Stress 




Delay (ns) 

Figure S7: Percentage of the total energy dissipation that is due to stress in the MMS mode as a 
function of switching delay. With increasing delay, we can apportion more of the switching role to 
stress and less to spin transfer torque so that the percentage of energy consumed by stress increases 
with delay The total energy dissipated of course decreases with increasing delay since spin transfer 
torque current is much more dissipative than stress. 
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Energy Dissipation with nickel and Terfenol-D 




Delay (ns) 

Figure S8: Energy dissipated in the MMS mode as a function of switching delay for Terfenol-D 
and nickel. The two curves are plotted for the maximum stress (40 MPa for Terfenol-D, and 100 
MPa for nickel) that can be generated in the magnetostrictive layer assuming that the maximum 
possible strain in the PZT layer is 500 ppm. Note that MMS mode becomes more energy efficient 
in Terfenol-D for delays longer than 6 ns. Thus, for delays exceeding 6 ns, Terfenol-D will be the 
preferred material; otherwise, nickel will be preferred between the two. 
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